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Page  4,  Eq.  (2.1)  should  read:  8q^/8n  =  q^/R^  • 

Eq.  (2.  5)  should  read:  a(s)  =  in  q^  +  in  [  R(s)  -  d(s)  ]  . 

Page  5,  13  lines  down:  "stramline"  should  be  changed  to  "streamline". 
Page  6,  3  lines  down:  Equation  number  (2. 12)  should  be  changed  to  (2. 11). 
Page  7,  8  lines  down:  Equation  number  (2.13)  should  be  changed  to  (2.12). 
Page  8,  Eq.  (3.12):  Superscript  hats  (^)  should  be  inserted  over  to  and  (i . 

Page  12,  4  lines  down  should  read:  "late  here,  in  a  slightly . " 

10  lines  down  should  read:  "bounded  and  7Q*  (cp,  f )  =  o(l/r). . . .  " 
Page  13,  2  lines  up  should  read;  "....  values,  and  Q  is". 

Page  14,  line  6  should  read: 

[  u,  v]=^^  ^  ^Vvdcpdt 

Page  14,  Eq.  (5.2)  should  read: 


4  [u,v  ]*  ^  ^  Vu7vdcpdi|( 

GV 

€ 

+  l/(2TTe^<nc)J  u  |  d  w  |  ^  v|dw| 


Page  15,  Footnote,  line  2  should  read: 
P 


f  u(cp)  dcp=  lim  (f  +  f  1 
J.ao  e  -  0  I  J-oc  J  e  J 


u(cp)  d  cp 


Page  16,  3  lines  up  should  read;  ". .  .  .  VAR,  using  (5.  6)  as  the.  . . .  ". 
Page  17,  1  line  down  should  read:  ".  .  . .  to  the  absence  of  any 
Page  21,  Eq.  (6.  5') should  read: 

f  tt/2 

A  [u  ]  =  D  [  u  ]  +  2  (a/iT )  \  u  (1,  2t)  dt 
o  Jo 

0 

Page  22,  Eq.  (6.  8)  should  read: 

f  n/2 

A  [  u  J  =  \  u  (1,  2t)  [  u  (1,  2t)  +  2a/  tt  3  dt 
J  °  °P 

0 

Page  36,  2  lines  down  should  read:  "....  numerically  using  the" 
Page  38,  2  lines  up  should  read:  ". .  .  .  property  (c)  it  follows.  ..." 
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SUMMARY 

The  injection  of  a  two-dimensional  jet  of  total  head  from  an  infinite 
plate  into  a  uniform  stream  of  lower  total  head  is  considered,  the  fluids  being 
assumed  inviscid  and  incompressible.  Steady,  irrotational  solutions  of  Euler's 
equations  are  found  for  Hj/H^  ^ao.  The  region  behind  the  jet  is  treated  as  a 
stagnant  wake  with  constant  pressure  equal  to  that  of  the  undisturbed  stream,  and 
the  jet  injection  angle  is  fixed  across  the  jet  opening. 

When  a  thin  jet  approximation  is  combined  with  Bernoulli's  principle,  a 
non-linear  boundary  condition  is  derived  along  the  vortex  sheet  separating  the  jet 
from  the  external  flow.  The  resulting  non-linear  potential  problem  (in  the  plane 
of  the  complex  velocity  potential)  for  the  external  flow  is  shown  by  Pal  to  be 
equivalent  to  a  variational  problem.  A  numerical  procedure  based  on  the  varia¬ 
tional  principle  and  the  Ritz-Galerkin  method  is  used  to  solve  for  the  case  of 
normal  injection  with  a  digital  computer.  The  pressure  distribution  along  the 
plate  upstream  and  streamline  diagrams  are  given.  ( 
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PART  I 
by 

Robert  C.  Ackerberg  and  Alexander  Pal 


1.  INTRODUCTION 


This  investigation  is  concerned  with  the  injection  of  a  jet  of  total  head  Hi  from 
an  infinite  plate  into  a  uniform  stream  of  lower  total  head  H2  (see  Fig.  la).  Two- 
dimensional,  steady,  irrotational  solutions  of  Euler's  equations  will  be  sought,  the  fluids 
being  assumed  inviscid  and  incompressible.  The  injection  angle  -a  is  fixed  across  the 
jet  opening,  and  at  large  distances  from  the  origin  all  motion  is  assumed  to  be  directed 
in  the  positive  x-direction.  The  region  to  leeward  of  the  jet  will  be  treated  a  a  stag¬ 
nant  wake  with  constant  pressure  equal  to  that  of  the  undisturbed  stream. 

Flows  of  this  type  occur  frequently  in  connection  with  VTOL  aircraft,  ground 
effect  machines  and  jet-flapped  wings.  The  basic  difficulty  in  their  analytical  study 
is  the  non-linear  boundary  condition  which  must  be  applied  along  a  streamline  sepa¬ 
rating  regions  of  different  total  heads.  To  ensure  the  continuity  of  the  static  pressure 
a  vortex  sheet  coincident  with  this  streamline  must  be  introduced.  The  extensive 
literature  dealing  with  this  subject  assumes  one  or  more  of  the  following:  1)  =  Ha , 

2)  the  jet  is  bounded  by  free  streamlines  and  walls  and  there  is  essentially  no  inter¬ 
action  with  the  external  flow,  and  3)  the  jet  injection  angle  is  small  and  linearized 
flow  is  assumed. 

Previous  work  using  this  model  for  90°  injection  has  been  carried  out  by 

12  3 

Taylor  ,  Ting,  Libby  and  Ruger  ,  and  Vizel  and  Mostinskii  .  Taylor  uses  a  rough 

theoretical  analysis  to  determine  the  shape  of  the  bounding  streamline  OC  near  the  jet 

exit  and  points  out  the  difficulty  of  verifying  any  results  experimentally  because  the  jet 

will  fill  a  wedge  of  nearly  40°  due  to  viscous  spreading.  Ting  et  al.  considers  the 

problem  when  Hj  is  slightly  larger  than  H2  and  uses  ordinary  perturbation  techniques 


to  correct  the  solution  obtained  for  equal  total  heads.  They  noted  the  difficulty  of  a 

straightforward  perturbation  procedure  due  to  a  non- uniformity  in  the  flow  deflection 

angle  along  the  streamline  OC  [near  the  jet  exit]  when  (H!  -  Ha)  -  0  .  Vizel  and 

Mostinskii  treat  the  problem  semi- empirically  by  introducing  a  coefficient  of  bending 

which  is  equivalent  to  specifying  an  average  radius  of  curvature  over  the  length  of  the 

jet.  They  compare  their  results  with  experimental  data  obtained  in  one  of  their  refer- 
t 

ences. 

The  fundamental  assumption  of  this  analysis  will  be  that  the  ratio  of  jet  thickness 
d  to  its  radius  of  curvature  R  will  be  small.  Physically  this  is  to  be  expected  when 
H1/Ha  -•  oo.  Using  the  simplified  flow  of  Fig.  1,  it  is  then  possible  to  take  the  basic 
non-linearity  of  the  problem  into  account.  * 

This  report  is  presented  in  three  parts.  In  Part  I  (by  R.  C.  Ackerberg  and 
A.  Pal)  a  non-linear  boundary  condition  is  derived  for  the  external  flow  based  on  the 
assumption  that  d/R  <  <  1  inside  the  jet  ($2).  Since  the  jet  shape  is  unknown,  the 
analysis  is  similar  to  potential  flows  with  free  streamlines  and  is  readily  carried  out 
in  the  plane  of  the  complex  velocity  potential.  In  §3,  a  non-linear  potential  problem  is 
derived  in  this  plane  for  the  logarithm  of  the  speed  along  any  streamline  in  the  external 
flow.  Formal  asymptotic  solutions  of  the  potential  problem  have  been  found  in  $4  which 
are  valid  close  to  and  far  from  the  jet  exit. 

In  Part  II  (by  A.  Pal),  the  equivalence  of  the  boundary  value  problem  with  a 

*  We  have  not  been  able  to  locate  a  copy  of  this  reference. 

$  4 

Motivated  by  this  study  Ackerberg  has  started  with  the  assumption  -*  oo 

instead  of  d/R  <  <  1.  The  two  assumptions  are  shown  to  be  equivalent  to  first 
order  for  all  streamlines  not  too  close  to  the  jet  exit.  See  the  footnote  on  p.  5. 


2 


variational  principle  is  discussed  in  In  '.'6,  a  numerical  method  of  minimizing 

the  variational  integral  based  on  the  Ritz-Galcrkin  method  is  given.  The  resulting 
system  of  non-linear  equations  is  solved  numerically  for  90°  injection  in  §  7  using 
Newton's  method. 

In  Part  III  (by  R.  C.  Ackerbcrg  and  A.  Pal),  the  numerical  results  are  compared 
with  experimental  and  theoretical  results  of  other  workers. 

Most  results  obtained  in  this  analysis  would  be  altered  significantly  by  viscous 
effects  arising  from  the  shear  layers  along  OC  and  AB.  In  addition,  this  flow  is 
notoriously  unstable  and  the  assumption  of  steady  flow  could  probably  not  be  maintained 
in  an  experiment.  Nevertheless,  it  is  expected  that  inviscid  results  will  be  valid  close 
to  the  jet  exit  where  a  determination  of  the  forces  on  the  plate  can  be  made. 


2.  THE  THIN  JET  APPROXIMATION 

Introduce  the  complex  space  coordinate  z  =  x  +  iy  with  origin  at  point  O 
(see  Fig.  la).  Since  the  locations  of  the  streamlines  OC  and  AB  are  not  known  in 
advance,  it  is  convenient  to  formulate  this  problem  using  a  coordinate  system  based  on 
the  streamlines  and  equipotential  lines.  Define  the  complex  velocity  potential  w  =  cp  +  i'K 
with  i)l  =  0,  and  i|/  =  m  >  0  along  OC  and  AB  respectively,  and  |  <  0  in  the  external  flow. 
The  vortex  sheet  along  OC  requires  that  two  velocity  potentials  be  used,  and  conditions 
in  the  jet  and  the  external  flow  will  be  denoted  by  the  subscripts  1  and  2,  respectively.  * 
The  flow  region  in  the  z-plane  maps  into  the  w-plane  as  shown  in  Fig.  2. 

Choose  a  curvilinear  coordinate  system  (s,  n)  in  the  jet  as  follows  (see  Fig.  lb): 

Since  4/  is  continuous  across  the  streamline  OC,  it  need  not  be  subscripted. 
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The  n-coordinate  of  a  point  P  in  the  jet  will  be  the  distance  PQ  measured  along  the 
cp  =  const,  line  through  P,  and  the  s-coordinate  will  be  the  arc-length  OQ  along  the 
bounding  streamline  OC.  The  condition  of  irrotationality  requires 

jj}  ^  A  ,  (2.1) 

c  r\  R  l 

where  Rx  is  the  radius  of  curvature  of  the  streamline  through  (s,n).  When  d/Rj  <  <  1  , 
we  may  with  good  approximation  assume 

Rj  ( s  ,  n  )  *  R  (  )  -  >1  ,  (2.  2) 


where  R(s)  is  the  radius  of  curvature  of  the  streamline  OC.  Physically  this  means  that 
the  streamlines  in  the  jet  crossing  an  equipotential  line  are  approximated  by  concentric 
circles.  Substituting  (2.2)  into  (2.  1)  and  integrating  yields 

£»  (v>-  t-  (1  (2'3) 

where  a(s)  is  an  arbitrary  function.  Along  the  free  streamline  AB  we  require 

n  )=  <1^.  ,  ft.  n  --(I  (<,)  .  (2.4) 


Substituting  (2.4)  in  (2.  3)  we  obtain 


Therefore, 


a  (  5  )  =  l' 


-/xr>  i 


i»Lr 


(  '»  ) 


(  •'> )  -  (U  v) 

R  ( *)  -  n 


(2.5) 


(2.6) 
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thus,  at  any  fixed  value  of  s  the  velocity  distribution  is  that  of  a  potential  vortex.  ^ 

The  conservation  of  mass  in  the  jet  requires  that  a  fixed  volumetric  flow  rate 
cross  every  section  s  =  constant,  i.  e. , 

d(‘ '  ) 

J  (|L(^h)cin--  (2>. 

c 

where  d^  =  d(oo)  is  the  uniform  jet  thickness  far  downstream.  Substituting  (2.6)  in 
(2.  7)  and  integrating  yields 


-[K(t;-ot(s>J  U  — <2‘8’ 


l.  e. 


cUO-d^lj  t 


(2.9) 


Thus,  the  thickness  of  the  jet  is  constant  to  first  order. 

Across  the  streamline  OC  the  pressure  must  be  continuous.  Using  Bernoulli's 
principle  this  requires 


h.V-**V*  -u'vhj  , 


(2.  10) 


where  the  A  denotes  the  velocities  at  some  point  along  the  bounding  stramline 

OC  [e.  g.  ,  qa  =  qx  (s,  o)]  .  Here  ?  are  the  total  heads  in  the  jet  and  the  external  flow, 

respectively.  When  (2.  10)  is  applied  at  an  infinite  distance  downstream,  a  relationship 


t  4 

When  a  t  n/2,  Ackerberg  has  shown  that  this  approximation  is  in  error  by  0(d/R) 
near  the  jet  exit  [|z|  =  O(d)]  as  a  result  of  neglecting  a  boundary  condition  along  OA. 
The  thin  jet  approximation  is  valid  in  the  region  |  z|  =  0(d/u),  where 
U  =  0J  Nevertheless,  the  motion  in  the  external  flow  (i.  e.  ,  the 

solution  of  Eqs.  (3.  1 3 ) - (3 _  17))  will  be  determined  correctly  to  o(l)  up  to  the  jet  exit 
because  the  inaccuracy  of  (2.  6)  modifies  the  boundary  condition  along  OC  only  to  a 
small  order. 
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between  the  basic  parameters  is  obtained.  Introducing  the  small  dimensionless 
parameter  u  =  (paq3  /Piqs  )  we  obtain 


H  i  ~  Ha 

i  n  ,  i 
Z 


(2.  12) 


Substituting  (2.  6)  with  n  =  0  into  (2.  10)  and  neglecting  terms  of  0[(d/R)s]  yields 


(2.  12) 


This  is  a  boundary  condition  along  OC  for  the  external  flow  which  relates  the  speed  and 
curvature  along  this  streamline. 


3.  THE  FLOW  IN  THE  EXTERNAL  STREAM1" 


In  the  external  flow  introduce  the  complex  velocity  potential 

I 

w  (z  )-  r  Lf  ,  (3' 

where  ♦  <  0,  -oo  <  C0  <  oo.  The  complex  velocity  is 

c { W'  •  -I  & 

■  (3-2) 

The  boundary  value  problem  for  the  external  flow  is  most  readily  formulated  in  the 
w  -  plane.  Introduce  the  logarithm  of  the  complex  velocity 

>-<*(*  y'  •  <3-3> 

Here  Q  and  9  are  conjugate  harmonic  functions  related  by  the  Cauchy-Riemann 


1  Hereafter  the  subscript  2  which  identifies  quantities  in  the  external  flow  will  be 
omitted  whenever  possible. 
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equations 


2a 

ot|/ 


'it 

o  G- 


(3.4a) 


(3.4b) 


and  defined  in  the  region  \|i  <  0,  -oo  <  cp  <  oo. 


The  curvature  of  any  streamline  is  given  by 


1 


c  e 


v  vy  . 

ROf. ~)  s  r  R 


da 


(3.5) 


where  we  have  used  (3.  4b)  and  the  relation  q  =  q  e  .  Applying  (3.  5)  to  the  bounding 

OOP 

r  «.  Q(cp,  0) , 

streamline  OC  and  substituting  in  (2.  13)  yields  [note  q  =  q  e  J 

^  003 


,  2otu*  ciiAJliQ  ) 

<-  1  u e  GfV-o 


Solving  for  dQ/di|i  we  obtain 

/  3(2  \  u 


m. 


0  1  << 


[c  e  a  Lf,c>j 


-—i  y*"  - 

^ 1 'if  ^  > 

Along  the  wall  DO  the  deflection  is  fixed.  Thus, 

0(f,  0;*0  Q&nij  Y^O^'^0 


(3.6) 


(3.7) 


(3.8) 


This  may  be  written  in  terms  of  Q  using  (3.  4b),  i.  e.  , 


9(j< 


—  (j  cl 


(.CM 


<j  f  -c ,  Y<c. 


(3.9) 


At  infinity  the  flow  must  be  uniform  and  undeflected.  Therefore, 
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0O  . 


(3.10) 


Q  (f,  y  )  0 ,  Q  (cf?  y  )  0  J  oh>  I  w  I  — ► 

Eqs.  (3.  7)-(3.  10)  are  satisfied  by  the  function  r(w)  =  0,  corresponding  to  a  uniform 
stream  and  for  which  the  flow  does  not  turn  through  an  angle  -a  near  point  O.  To 
satisfy  this  condition  F  must  have  a  logarithmic  singularity  at  n  =  0  (i.  e. ,  w  =  0) 
corresponding  to  a  stagnation  point  flow  in  a  corner.  Thus 


r(w)-  i  U(en^)t  r*(v)  ,  (3-n) 

where  0  >  arg  w  >  - rr,  and  T  (w)  is  bounded  at  w  =  0. 

On  introducing  the  non-dimensional,  scaled  complex  velocity  potential 

'V  g  c friof-  -  2  (3.12) 

s  1  t  ^  * 

in  Eqs.  (3.  7)- (3.  11),  the  following  non-linear  potential  problem  is  obtained  for  Q2 : 


C 

— t— > 

A 

Vt 

•A 

-  *3  <  <  -XJ  / 

(3.13) 

;-(x 

A 

0  for  y  -  O 

*\ 

O 

V 

br- 

(3.  13) 

c )Q 

af  " 

—  in  Q  p  Or 

f --o^  y>q 

(3.15) 

Q(<f 

VV  | 

(3.16) 

bounded  at  t  =  0,  and 


4.  ASYMPTOTIC  EXPANSIONS  OF  T(w) 


Formal  asymptotic  expansions  of  T(w)  which  are  analytic  and  which  satisfy 
the  boundary  conditions  term  by  term  without  inconsistency,  were  constructed  for 
w-0,  and  w-  oo,  Denoting  e^1  w  =  W ,  we  obtain  for  W  0, 

ri-W)~a&x\vff  (4.1, 

n "0  m 

where  d  isingmeial  a  polynomial  in  fnW  with  real  coefficients.  Coefficients  have 

been  obtained  recursively  up  to  n  =  2  with  no  inconsistencies  provided  a  single 

logarithmic  term  W2fnW  is  introduced.  The  coefficients  which  multiply  integral 

powers  of  W  are  indeterminate  by  this  formal  procedure. 

If  a  =  n/2,  the  following  simpler  expression  is  obtained1 
A  °"C 

r(-w)~  i  W+  1  anW  <4-2> 

rt-o  } 

'aO 

where  a0  and  a2  are  undetermined  and  a:  =  -e  .  From  the  numerical  computations 

inPartIIfora  =  tt/2,  aQ  =  0.0450...,  ax  =  -0.955...,  and  ao  =  -  1 .  58.  .  .  Note 

-ao 

that  the  relationship  ax  =  -e  is  satisfied  with  good  approximation. 


In  the  case  W  -  oo,  the  form  of  the  formal  asymptotic  expansion  does  not 


depend  on  a.  We  find 
A 


OO 


r(-w)~£  p^^w)\\r'ntl/±> 


(4.  3) 


hzO 


9 


mm 


i 


;« 


f 

> 

3 


where  Pn(L)  denotes  a  polynomial  (which  seems  to  be  of  o-der  n)  with  real  coef¬ 
ficients.  From  the  first  few  terms  it  appears  that  the  P  's  can  be  obtained  recur¬ 
sively  with  the  exception  of  their  constant  terms.  Thus, 

r(-W)~  A  c  VV  +  (  "iT  A,  fn  H/  '■A,)VV  (4.4) 


* 

Numerical  computations  for  a  =  tt/2  yield  A  =  -0.  797.  .  . 

The  equations  of  the  streamlines  may  be  found  by  integrating  (3.  3),  i.  e.  , 


-fc. 

cl 


2  -- 


(4.  5) 


Putting  w  =  <p  >0  in  (4.  2)  and  (4.4)  and  substituting  these  results  in  (4.  5),  we  obtain 
the  asymptotic  behavior  of  the  streamline  OC  close  to  and  far  from  the  jet  exit.  These 
results  are  readily  expressed  in  terms  of  the  nondimensional,  scaled  complex  space 
coordinate 


z- 


>  v 

L  I 


(4.6) 


where  z  denotes  the  complex  coordinate  of  a  point  on  OC.  For  |z|  0  in  the  special 

case  a  =  tt/2  t 


Although  these  expansions  are  derived  in  a  purely  formal  manner,  theoretical 
results  of  Pal^  suggest  that  they  are  indeed  asymptotic  to  T(w).  It  is  proved 
there  that  Q  is  asymptotic  to  the  leading  terms  of  the  expansions,  and  formulas 
are  obtained  expressing  d0o  and  Ao  with  integrals  involving  the  boundary  values 
of  Q.  Note  that  Q(cp,  0)  =  0(cp'3/? )  for  cp  —  +  oo  whereas  T(w)  =  Odwl*1/2)  for 
| w  |  -  ro  in  agreement  with  the  formal  asymptotic  expression  (4.  4). 
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(4.  7) 


x-ifs  o(vA), 


and  for  z  "*  a>, 


Using  only  the  first  term  of  (4.  7),  and  the  first  term  of  (4.  8)  with  an  experimentally 
determined  constant  added,  these  asymptotic  formulae  are  plotted  along  with  the 
numerical  result  for  a  =  n/2  in  Fig.  3.  The  asymptotic  curves  fare  in  well  with  the 
numerical  result. 
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Alexander  Pal 


The  author  is  indebted  to  Dr.  R.  C.  Ackerberg  for  his  cooperation  in  the 
research  project  out  of  which  this  section  grew.  Frequent  discussions  with 
him  of  details  of  this  work  and  his  careful  review  of  this  section  are  con¬ 
sidered  especially  helpful.  He  also  helped  by  programming  some  of  the 
subroutines  used  in  the  numerical  work. 


5.  VARIATIONAL  PRINCIPLE 


In  Part  I  of  this  paper  it  was  shown  that  the  fluid  flow  problem  of  the  interaction 
of  a  fast  jet  with  a  potential  flow,  which  is  a  parallel  flow  at  large  distance  from  the 
orifice  can  be  reduced  to  the  solution  of  a  plane  boundary  value  problem  which  we  formu¬ 
late  there,  in  a  slightly  different  way,  using  reflection  on  the  cp-axis. 

Let  G  denote  the  whole  (cp,  ♦  )  plane,  slit  along  the  -cp-axis.  Then  a  function 
Q(cp,  \|f )  has  to  be  found  which  satisfies  the  following  conditions: 

Conditions  in  G 

(a)  Q(cp,  t|[ )  is  harmonic  in  G  and  symmetric  to  the  real  axis. 

$  i 

(b)  The  function  Q  (cp  ,  ^  )  =  Q(cp,  |  )  -  (a/n)  £n  r  (r  =  (cp2  +  \Jf3  ) 8  )  is 

bounded  and  VQ  (cp,  ^  )  =  o(r)  uniformly  as  r  -  0. 

(c)  Q(cp,  ilf )  -•  0,  V Q(cp,  \|f )  =  o(l  /  r ) uniformly,  as  r  -  oo. 


Boundary  Conditions 

(d)  Q  and  are  continuous  on  the  +cp-axis.  Here  n  is  the  outer 

o  n 

normal  to  the  domain  G  on  f  =  0,  i.  e.  ,  ^  ^  on  the  lower  side  and  =  -  ^  on  the 

on  A  V  Of 

upper  side  of  the  cut. 


(e)  ^  =  -  sinh  Q  on  ♦  =  0  ,  ^  >  0,  (5.1) 

on 

(f)  Q(cp,  0)  =  (%  3^?)  as  ^p-  +oo. 


Conditions  (b),  (c),  and  (f)  are  more  stringent  than  necessary.  It  was  shown 
that  the  problem  has  a  unique  solution  under  much  milder  conditions  (see  also  footnote 
on  page  16).  However,  since  this  solution  satisfies  the  stricter  conditions  above  (including 
(b),  (c),  and  (f)  ),  we  may  use  this  formulation  for  convenience.  It  should  also  be  noted 
that  the  formal  asymptotic  expansions  obtained  in  Part  I  are  consistent  with  these  condi¬ 
tions. 
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Functions  which  satisfy  conditions  (b),(c),(d),  and  (f),  (but  are  not  necessarily 

harmonic  or  satisfy  the  boundary  conditions)  will  be  called  "admissible",  and  the  family 

of  admissible  functions  with  a  fixed  value  of  a  will  be  denoted  by  Q  . 

a 

An  important  property  of  the  family  is:  If  the  functions  u,  v  belong  to  Q,  , 
then  for  all  X(0  <  X  <  1), 

vv  =  A  +■  ( i "  A  ^  v 

also  belongs  to  G  This  maybe  expressed  concisely  as  follows:  The  function  space 
U  is  convex.  The  statement  is  obvious;  to  see  it  one  only  needs  to  verify  that  w  satis' 
fies  all  criteria  of  admissibility  if  u  and  v  do. 

The  main  difficulty  in  solving  the  above  boundary  value  problem  stems,  of  course, 

from  the  non-linearity  of  the  boundary  condition  (5.  1).  The  fact  that  Q  is  harmonic  sug- 

6*# 

gests  a  variational  method  based  on  the  Dirichlet  principle. 

Generalized  Dirichlet  Integral.  We  first  extend  the  notion  of  the  Dirichlet  integral 
slightly  to  admit  functions  which  have  a  logarithmic  singularity  at  the  origin,  such  that 


The  name  convex  is  suggested  by  the  analogy  to  finite  dimensional  ««Daces;  e.  g.  ,  a  ball 
is  convex  because,  if  the  points  p,  q  belong  to  the  ball,  then  the  entire  straight  segment 
connecting  p  and  q  belongs  to  the  same  ball. 

The  Dirichlet  principle  states  that  in  a  domain  G  with  a  sufficientlv  smooth  boundary, 
there  is  a  function  Q  which  makes  the  Dirichlet  integral 

D[Uj  -■  i  Jf(VQ)ictyoty 

Cr 

minimum  under  the  constraint  that  Q  assumes  given  continuous  boundary  values,  Q  is 
harmonic  in  G. 
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they  satisfy  condition  (b). 


domain 


We 


will  use  the  following  notations: 

1.  will  denote  the  subdomain  e  <  r  <  X  of  G; 

c 

2.  [u]  ,  [u,  v]  will  denote  the  Dirichlet  integral  of  u  over  the 

€  C 

,  and  the  associated  bilinear  form,  i.  e. 


0/  l«  J  ;  |-y  [U  ,  u  J 


3.  0  subscripts  and  oo  superscripts  will  be  omitted  in  these  forms,  as 

from  A*  defined  below, 

e 

4.  C ^  will  denote  the  circle  r  =  X  . 


The  Dirichlet  integral  is  for  admissible  functions  in  general  divergent.  In  order 
to  separate  its  "finite  part",  we  define 


and 


2  *  [<<J  s  ,a /[_<.< ,  i<  J  • 

It  is  easy  to  see  that  if  u(r,  y),  v(cp,  i|i  )  satisfy  conditions  (b)  and  (c),  then  A^  Lu,  v]  has 
limit  a[u,  v  J  when  €  -  0  and  X-oo,  (and  similarly  ‘^[uJ-aLu]  =  A[u,  u]). 

If  all  integrations  in  the  preceeding  definitions  are  extended  only  to  the  upper 

(lower'  halfplane,  we  can  define  the  functionals  ,  Lu],  A  ,  [u  ]  ,  etc.  (A^  [ u 3 ,  A  [u] 

*  >  '  '  £|  ~ 

etc.  )  completely  analogously.  Thus 


We  mention  some  important  properties  of  the  functional  aLu]  : 

(1)  Let  the  functions  u(cp,  iji),  v(cp,  t  )  be  admissible,  say,  u  c  d  ,  v  e  Qn . 

a  d 

Then  u  +  v  e  (l  .  Q  and 
a  +  p 

A  [  u  f  v  J  '  A  l.“  J  '  A[  o. ,  yj  ]  f  A  [  V  J  (5.  3) 

Eq.  (5.  3)  is  an  immediate  consequence  of  the  identity 

Z.c  |_U  t-  v  J  -  Ac-  l  U  J  f  J  /\  l  u  ,  s/  J  t  Ac  W  J  , 

if  we  observe  that  A  [u],A  [v],  A  Lu  +  v]  all  have  limits  as  €-0. 

£  £  C 

(2)  Suppose  h(cp,  t  )  is  a  harmonic  function  in  Q,^ ,  and  k(cp,  ^  )  is  a  function 
in  (thus  k  has  a  finite  Dirichlet  integral  in  G).  Then, 

*v,kJ-  f’i  h  (-f,C  K  .  *  (5-4) 

-  A' 

The  proof  of  Eq.  (5.  4)  can  be  found  in  Appendix  1. 

Substituting  (5.4)  into  (5.  3),  we  obtain 

oo 

A[h^]-Ath. j  p|  h^r,c)k(f,0)d^i  D|>J  .(5.5) 

-  OO 

(3)  The  functional  A[u]  can  be  defined  alternatively  by  ordinary  Dirichlet 
integrals.  Let  X  be  any  positive  number.  Then 

Alu  J  '  Dj  <4  J  T  D*L  u*  J  |  Ct</S  -  JJJ.  Ul)  (5.6) 

C 

where  u  =  u  -  (a /r )  / n  r.  The  proof  of  this  identity  can  be  found  in  Appendix  2. 

Formulation  of  Variational  Problem.  If  u  is  admissible,  then  the  boundary 

integral 


P  denotes  here  the  "Cauchy  principal  value"  of  the  integral: 

V  J  a  ( tf  ) ot  Vf  =  Ztm  [  C  *  *  f  l 

t  >  t-  J  OC  J  c  J 
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(5.7) 


8  j_  a  ]  j  [  Ccil\  *  i  (  y,  i< )  -  1  j  ^  '•(' 

o 

is  always  convergent.  From  A[u]  and  B[u]  we  form  the  composite  expression 

F[  a  J  a]  f- 6[“J  .  (5.8) 

We  mention  here  a  property  of  F[u]  which  will  be  useful  later:  F[u]  is  a 

"convex"  functional;  i.e.,  if  u  e  (l  ,  v  c  (1  ,  and  for  some  \  ,  0  <  X <  1 ,  w  =  \u+(l-\)v, 
-  a  a 

u^v,  then 

^  [  vvj  <  AFl  *■(>-))  Flv]*  (5.9) 

The  proof  of  this  property  can  be  found  in  Appendix  3. 

We  consider  now  the  following  variational  problem  (problem  "VAR"):  Find  an 
'admissible'  function  Q(cpf  <|i)ef2a>  such  that  fLQ]<  F[Qj]  for  any  e  Q  different  from 
Q.  Although  the  existence  of  a  function  Q  will  not  be  established  in  this  paper,  we  will 
show  that  the  problem  VAR  is  equivalent  to  the  problem  BV;  thus  if  problem  BV  has  a 
solution,  it  is  also  a  solution  of  problem  VAR,  and  conversely.  Such  a  solution  is  by  its 

■ff. 

minimum  property  necessarily  unique. 


Equation  (5.  9)  is  obviously  satisfied  by  an  (ordinary)  function  F(u)  which  is  convex 
from  below.  Moreover,  convex  functions  can  be  defined  by  (5.  9). 


r*  j 

The  existence  of  a  solution  of  the  problem  VAR  is  shown  (see  Pal  )  in  a  stronger 


form.  It  is  proven  there  that  there  is  a  unique  function  Q  admissible  in 
sense,  which  renders  F[Q]  the  smallest  in  a  much  wider  function  class 

requires  for  membership  D  [u  ]  <  oo  ,  D^[u]  <  oo  and  B[u]  <  oo,  instead 
properties  (b),(c),(d),  and  (f). 


the  above 

V  <aa 


A  different  proof  of  the  equivalence  of  BV  and  VAR,  using  (1.6)  as  the  defining  equa¬ 
tion  of  a[u]  was  given  by  Ackerberg.  ^ 


I  am  thankful  for  this  observation  on  the  uniqueness  to  R.  C. 


Ackerberg. 
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We  call  attention  in  this  variational  problem  to  the  absense  of  any  kind  of  boundary 
conditions  along  the  real  axis.  The  fulfillment  of  the  boundary  condition  will  be  shown  to 
be  a  consequence  of  the  minimum  property  of  the  solution.  Similar  situations  are  not 
infrequent  in  variational  problems.  ("Free  boundary  values",  "natural  boundary  condi¬ 
tions")^. 


Proof  of  the  equivalence  of  the  problems  BY  and  VAR 

I.  We  first  show  that  if  the  function  Q(cp,  i|f )  satisfies  the  conditions  of  VAR,  then 
it  is  also  a  solution  of  BV. 

(a)  The  function  Q  must  be  harmonic  in  G.  If  it  is  not,  then  there  is  a 

circular  disk  uu  entirely  in  G,  in  which  Q  is  not  harmonic.  Then  it  is  possible  by 
Poisson's  integral  to  construct  a  function  such  that  Q=Q^  on  the  boundary  and  out¬ 

side  uu ,  and  Qj  harmonic  in  x.  Then  by  Dirichlet's  principle, 

DJQJ<  DJGU  . 

(D^  denotes  here  the  Dirichlet  -  integral  in  x).  If  c  is  chosen  so  small  that  Cg  is 
completely  outside  x, 

[G  J  <At  CQ] 

since  Q=Q^  outside  x.  Hence 

A[C.(]  <  ^tCv]  , 

a  contradiction  to  the  assumption  that  Q  satisfies  the  conditions  of  problem  VAR. 

(b)  The  function  Q  must  be  symmetric  to  the  cp-axis.  Suppose  e.  g.  that 


A,  La]  a  a.[g]  . 


Then  we  define  the  new  admissible  function 


QL(<f,v ) r 


Q(f ,  -y ) 


>f  y  >  C  ■  J 
<f  V  <  C 
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Clearly 


A  L U1  ]  ^ 

which  by  assumption  is  only  possible  if  Q^  =  Q. 

(c)  We  will  use  the  definition  of  the  first  variation  of  functionals  customary 
in  variational  calculus.  If  6  Q(;p,  i(f  )  is  an  arbitrary  variation  of  Q,  and  $[Q]  is  a 
functional,  then 

£  tc,  ]  --  |  ^  ‘f  l  *  1  A  oU  J  _  c  (5. 10) 

is  the  "first  variation"  of  $[Q].  Thus  we  show  that 

£FLQ]-  pj  [Q  (vf ,  c )  +  v(if )  im/x  Q  C^cjJ  5 

^  (5.H) 

where  T](cp )  =  1 ,  if  cp  >  0  and  r](cp)  =0  if  cp  <  0. 


Equation  (5.  11)  immediately  implies  that  if  Q  is  a  solution  of  the  problem  VAR, 


then 


(Lw  )- 


c 


if  V  -  I ' 


G(r,0  'f  <f  >c  , 

i.e.  ,  Q  satisfies  all  conditions  of  BV. 

Equation  ( 5.  1 1 )  is  an  easy  consequence  of  ( 5.  5).  Substituting  h  =  Q,  k  =  \6Q  in 
(5.  5)  and  taking  into  account  the  symmetry  of  Q,  6  Q, 

X.-1 

a  [a  c  A6'&]  -  a  f  <f.t  )<$■(;  (v,c;+  /nisei]. 

O 

This  implies  by  definition  (1.  10)  that 

/,K_’ 

Sa[Q]-  (  Gt (f,c; .  (5.i2) 

'o 

On  the  other  hand,  from  the  definition  of  bLq]  (Eq.  (5.  7)  ), 

o  6  [Q  ]  -  j  Uv.k  G  (  c )  o  G  ( ^ ,  l  ,> 


However,  in  the  present  context,  only  such  variations  6Q  are  admissible  for  which 

Q  +  \  6Q  c  0  for  all  sufficiently  small  values  of  \|.  For  this  it  is  sufficient  if 
a 

bQed^.  Furthermore  if  is  sufficient  to  consider  only  variations  symmetric  to  the  cp-axis. 
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is  obvious.  Combining  with  (5.  12),  we  obtain  (5.  11). 


II.  Suppose  now  that  Q(cp,  i|»  )  is  a  solution  of  problem  BV.  Then  by  (5.  11), 
6F[Q]=0f  i.  e.  ,  F[Q]  is  stationary.  To  show  that  F[Q]  has  an  absolute  minimum  in 
Q,  we  will  use  the  "convexity"  of  the  function  space  and  of  the  functional  F. 

Let  Qj  i  Q  be  any  other  admissible  function.  Consider  the  functions 

■5  (  f  'V  ’  U  ■=•  ( /•  A  )  Q  *  >  (vi  •  (  c  <  b  S  1  ) 

By  the  convexity  of  Cl  ,  S  also  belongs  to  Cl  .  Denote  F[Sl  =  i{\).  6  F[Q]  =  0  implies 

a  a 

'(())=  0.  However,  by  the  convexity  of  F[u],  the  function  H \)  is  convex  from  below 
in  (0,  1);  thus, 

FLuj  ]*<HD  >f(C.V=-  f  LG)  , 

i.  e.  ,  F[u]  has  an  absolute  minimum  when  u  =  Q. 
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6.  CONFORMAL  MAPPING  INTO  A  CIRCLE 


In  the  numerical  solution  of  the  variational  problem  VAR  it  is  desirable  that  the 
domains  of  integration  be  finite.  Therefore,  the  half-plane  Y  <  0  is  conformally  mapped 
into  the  half-circle  disk  G  •  s  :  |s|<l,  Ims>o|  of  the  complex  s-plane.  The  mapping 


function  is 


and  its  inverse 


r  /  \  i  -  w  )  ^  f-  i 


(6.  la) 


(6.  lb) 


Here  (-wp  is  defined  by  providing  the  w-plane  with  a  slot  along  the  positive  imaginary 
axis;  thus,  0  <  arg  (-w)^  <  . 


( see  Fig.  4  ) 


iT 

Let  s=pe  .  The  mapping  of  the  positive  cp-axis  into  the  half-circle  c=l,  0  <  t  <  rr  is 
described  by 

Cf  ~  [cot  {  r /■?-  >  *  (6.  2) 


20 


Evaluation  of  F[u]  in  the  s-plane.  To  take  the  logarithmic  singularity  of  the 
function  T(w)  =  Q  -  i  P  at  the  origin  of  the  w-plane  into  account,  we  introduce  the  functions 

R(S)  -  -J*  .[\\  --f-'  '  (6.3) 

j  l 

i' ;  t  2  ('  f  <r'»7  i  I 

f<x  Rb)-  ?  ^  4.  '  (6.3') 

Write  I'(s)  =  F(W(s))  and 

r(s>- R(s>  ,  ,6.4, 

a  H,r)-=  C«,r)-  +(s,,r  J  •  ,6.4') 

G 

In  general,  for  any  admissible  u(o»T  ).  the  function  u^(p, t)  =  u  -  r(p,  t)  is  bounded  in 

p  <  1.  Further  conditions  on  u  are  specified  below. 

—  o 

The  functional  A  [u]  can  be  transformed  into  the  s-plane  using  (5.  5).  Replacing 

there  h  by  r(p,  t  ) ,  k  by  u  and  observing  that  r  (1,  t)  =  ql/tt,  we  obtain 

°  IT 

4|>J  =  D[uJ  ,  2  j  )dr  -  (6.5) 

0 

Since  the  Dirichlet  integral  is  invariant  with  respect  to  conformal  mapping,  the  term 

D[Uq]  may  be  evaluated  in  the  s-plane.  A[r],  a  constant,  is  irrelevant  with  respect  to 

the  variational  problem.  Therefore  we  will  use  the  form 

•/z. 


where 


[uj  -  Dl“c 

]  >  )  uc(i,jt )eft  , 

(6.5') 

* 

tl 

N 
* — 1 
£ 

J  *loj. 

(6.6) 

:r/j 

>  r  [  c^u  ** 

. ) 

'zrq  dt 

0 

ejo[  i  uc  (  1 9 

)J  US  t  -  «  p[~A  u°( 1  > 

>  Wr 

'  (6.7) 
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When  uq  is  harmonic,  (6.  5'  )  can  be  expressed  in  the  form  of  boundary  integrals  by 
transforming  D[uq]  with  Green's  formula: 

rr/St 

ALUJ  -  j  ue(l,2t)[  +  l]  dt  .  (6.8) 

o 

Boundary  conditions  at  s  =  + 1  and  s  =  -  1 .  Q(cp,  t)  is  expected  to  possess  properties 
(b),(c)  and  (f)  of  "admissible"  functions.  In  the  s-plane  this  can  be  expressed  as  follows, 
(b')  Qq(1  ,  rr)  =  0,  and  7Qq  =  o(|s  +  1  I'1) 

as  s  “•  -1.  (Here  as  above.  Q  (p,  t)  =  Q(P,  t)  -  r(p,  t).  ) 

o 

(c  *)  Q  (1,0)  =  0 ,  and  VQ  =  o  ( |  s  -  1  |  * )  as  s  +1. 
o  o 

(f7)  Q(1,t)  =  0(t3)  (6.9) 

as  t  “*  0.  This  implies 

qo(1.t)  =  ^t2  +  °(t3)  (6.10) 

as  t  “*  0. 
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7.  APPLICATION  OF  THE  RITZ-GALERKIN  METHOD 


The  variational  problem  VAR  was  solved  in  the  special  case  a=n/2  (normal 
injection)  numerically,  with  the  Ritz-Galerkin  method.  This  method  is  often  used  with 
advantage  when  the  value  of  the  variational  integral  is  the  most  important  result,  e.  g.  , 
in  problems  involving  eigenvalues  (elastic  buckling,  resonant  frequencies  in  vibration 
problems).  It  is  known  from  experience  that  usually  the  value  of  the  variational  integral 
can  be  determined  much  more  precisely  than  the  minimizing  function.  In  the  present 
case  only  the  minimizing  function  is  of  physical  significance.  Nevertheless,  the  Ritz- 
Galerkin  method  (referred  to  by  R  -  G)  was  found  quite  suitable  and  probably  better  than 
a  relaxation  type  method.  Our  choice  fell  on  the  former  for  the  following  reasons: 

(1)  In  the  R-G  method  the  unknown  function  is  approximated  with  a 
linear  combination  of  known  functions  (we  will  say  that  these  functions  form  a  "base"). 
The  base  can  be  chosen  in  such  a  manner  that  they  contain  much  of  the  information  known 
about  the  solution  before  the  computation  is  started.  In  the  present  problem  it  is  known 
that  the  solution  is  harmonic,  symmetric  to  the  real  axis,  and  smooth  on  the  boundary. 
The  base  can  be  chosen  to  satisfy  these  conditions  identically.  Boundary  conditions  at 
s  =  - 1  and  s  =  +1  can  be  taken  into  account  by  proper  combination  of  the  elements  of  the 
base.  Thus  a  significant  portion  of  the  computing  work  can  be  saved  if  known  properties 
of  the  solution  are  taken  into  account  by  a  judicious  choice  of  the  base. 

In  the  present  problem  the  R-G  method  requires  the  determination  of  a  one- 
dimensional  array  of  coefficients,  (p^  in  Eq.  (7.  2))  as  opposed  to  the  two-dimensional 
array  of  grid  points  in  the  relaxation  methods.  (Actual  computation,  discussed  below, 
showed  that  only  15  coefficients  yielded  as  satisfactory  results  as  probably  thousands  of 
values  in  grid  points  would  have  done.  ) 
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(2)  The  circular  shape  of  the  domain  makes  the  use  of  rectangular  grids 
clumsy.  On  the  other  hand,  if  the  w-plane  is  chosen,  the  infiniteness  of  the  domain  is  a 
source  of  trouble. 


According  to  the  R-G  method,  we  use  a  formal  expansion 


r„u)  =  res) -  u  ~- 
G„u>  i  u 


'  f  h. 

*  0 

^  f  fell  f  / 

~~  4 


(7.1) 


~  Z  ,  (7.2) 

where  H  (s)  are  analytic  functions  to  be  specified,  h  (n,  t  )  =  Re  H  (s),  and  the  p  are 
n  ’  r  n  n  rn 

unknown  coefficients.  To  take  full  advantage  of  the  economy  of  the  R-G  method,  the 
functions  hn(p»  T)  will  t>e  chosen  such  that  they  satisfy  the  following  conditions: 

1)  h  (p,  t  )  is  harmonic  in  p  <  1,  continuous  on  p  =  1. 

n 

2)  h  (p,  t  )  is  symmetric  to  the  real  s-axis. 

n 

3)  hn(  1 ,  T  )  must  form  a  complete  system  in  (0,  n),  i.  e.  ,  every  square 


integrable  function  u(t)  can  be  expanded  in  a  series 


u(T)r  Z 

A 


h  ( 1 ,  t )  need  not  be  differentiable  in  T=0,rr,  since  Q  (1,t)  is  not  differentiable  there, 
n  o 

The  functions  h  (p,  t)  will  be  determined  such  that  on  the  half-circle  p  =  1, 
n 

0  <  t  <  tt 


This  expansion  need  be  valid  only  with  the  exception  of  a  set  of  measure  zero. 
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n  (  / 


(  r>  b  n  i 


{n-  Gy  i  . .  )  , 


(7.  3') 


h 2n-l(  '  '  1  )  = 


it  uni 


(  k\  I  >  si.  /  .  i  >  )  • 


(7.3") 


Q^(l,  T)  can  then  be  expanded  in  an  infinity  of  ways  into  the  series 

Cv0  (  h  T  )  ~  ^  (/>„*  CcMil  t  p''  itvin?) 

which  is  valid  in  0  <  T  <  n. 


(7.4) 


Conditions  1)  and  2),  combined  with  (7.  }>')  and  (7.  3"),  are  sufficient  to  determine 


h  (p  ,  T  )  uniquely, 
n 

,'t  )=  tf*  trinl  (ti'tjIjA,...)  (7.6) 

is  obtained  quite  straightforwardly.  The  functions  h1  ,  are  somewhat  more  cumber- 

Zn  -  1 

some  to  compute.  Taking  into  account  property  2), 

{!<»>  tit  if  C  <  A  , 

r 

-  J.n  n(  if  7T  <  r  ^  2w  , 

By  Fourier  series  expansion  of  h  j(1,t),  and  matching  with  harmonic  functions 


term-by-term, 


*XJ  \  K 

is  >  t  —  ^  > 


(7.7) 


(n=l,2,...),  where  ?.  =  -r,  =  1;\  =0  forevenvaluesof  v  and  \  =  1  for  odd 

o  2  12  v  v 

values  of  v. 


Actually,  a  function  $(f)  with  continuous  §  ' (T  )  in  (0,  n),  can  be  expanded  into  series  of 
the  forms 


$  ("i  )  "  2.*  (\  „  fel  C 

f  (r;  -  2-T  ^  UK  nT  ■ 


(7.  5  ') 
(7.5") 


The  redundancy  (and  consequent  greater  flexibility)  in  the  expansion  (7.  4)  allows  a  closer 
approximation  of  Q::(  1 ,  T  )  then  would  be  possible  with  series  of  either  the  type  (7.  5  /  )  or 
of  (7.  5").  The  latter  series  converge  in  the  neighborhood  of  the  endpoints  0,  n  rather 
slowly  in  general,  because  either  all  the  base  functions  or  all  their  first  derivatives  van¬ 
ish  at  the  end  points.  Therefore  the  approximation  of  a  function  with  non-zero  values 
and  fi r st  de  ri vatives  of  0  and  r  is  not  practical  with  the  series  (7.  5  ')  or  (7.  5"). 
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The  corresponding  complex  forms  are 


- 


n 

S 


y 

K 


oo 


€ 


A  h-K 


U  -  o 


r\‘ 


j. 


<  *>  =  0,  1,2,...  ) 
(  n  *  I  ,  2  j  .  t  ) 


(7.8) 

(7.9) 


The  expansions  (7.  7),  (7.  9)  converge  too  slowly  if  p  is  close  to  unity.  Therefore,  it 
was  necessary  to  use  the  following  alternative  expression: 


1  {ft**-  s'")*.  77T  +  Z±  ( ^  J 


(7.10) 


«  =  / 


Equation  (7.  10)  was  obtained  from  (7.  9)  by  a  tedious  rearrangement  of  the  series.  The 

equivalence  of  (7.  9)  and  (7.  10)  is  simpler  to  show  by  substituting  the  power-series  expan- 
1  ^  s 

sion  of  -tn  — —  into  (7.  10).  The  two  series  are  complementary  i.  e.  ,  (7.  10)  is  suitable 
for  computation  when  |  s  |  is  not  too  small,  and  (7.  9)  should  be  used  for  small  |  s  | . 


Variational  integrals.  Approximating  Qq  with  a  finite  number  of  terms  of  the 
series  (7.2),  both  A[Q]  and  B[Q]  become  functions  of  the  coefficients  p  .  We  write 


and 


J  =  <•"(  ) , 

P  (  P  v  /  ■  '  ■  )  P  r. i-|) 

F[G]=  S(p)r{i{p) 


(7.  11) 


The  functions  6,8  can  be  obtained  from  (6.  3  )  and  (6.  5)  respectively,  by  the 


substitution  of 


(7.  12) 


n  =  0 


N  -/ 


Thus,  6( p)  is  a  mixed  quadratic  form, 

S(P )={Z  t  bnp 


rn-0  *~c 


(7.13) 


The  coefficients  a 


,  b  are  defined  by  the  formulas 
mn  n 


a 


*»i  a 


S  A*(  /,t  )Kh (t  )dr  , 

bti  -I  \  A„(  1,1  )dr  , 

c 

and  K  (T)  =  ["  -  h  (p,  t)  ,  .  From  Eqs.  (7.  8)  and  (7.  10) 

n  L  o  P  n  J  _  , 

p  =  1 


K 


K 

)  -  n  t’M  n  T 

> 

2  m 

)  -  i  . 

f  h  cts  m  T  hi  - 

t 

n 

l  2- 

i  1 

+  2  l 

X 

H  M  T 
il  n  l 


k  I 


(7.14) 

(7.15) 


(7.16') 


(7.  16") 


The  matrix  elements  a  ,  b  can  be  expressed  from  (7.  14)  and  (7.  15)  in  closed  form. 

mn  n  r 

Since  this  involves  a  considerable  amount  of  transformation  of  series,  and  is  quite  labo¬ 
rious,  the  details  are  omitted.  The  closed  formulas  obtained  for  the  matrix  elements  are 
very  useful  in  the  numerical  work,  and  are  listed  in  Appendix  4.  These  formulas  are 
consistent  with  the  symmetry  of  the  matrix  a^^,  and  were  checked  by  numerical  integra¬ 
tion. 


No  similar  reduction  of  3(p)  is  possible  because  of  the  non-linearity  of  the  boundary 

conditions.  Thus,  from  (6.  5) 

„  JT/_i  .  , 

fh  (p  )  *  I  E(tiP')  ''[Cc\t.  L((  ,p)  -  *t  cit  (7.17) 

where  i  % 

E(i,p)  =  t<P  (r  V  •  <7‘18> 

The  problem  of  finding  the  function  Q(.~,t),  for  which  FLQ]  assumes  its  minimum  is 
thus  reduced  (in  approximation)  to  that  of  finding  the  minimum  of  the  function  of  N  vari¬ 
ables  $(p)  =  6(p)  +  6(p). 
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The  function  4>(p)  will  be  minimized  by  constructing  a  sequence  of  "admissible" 
coefficient  families  |p^j(r  =  1,2,3,...)  such  that 

!nf$(p)  »■ 

In  the  numerical  work  for  the  determination  of  the  function  Q,  the  coefficients  p^were 
subjected  to  the  constraints  given  below.  A  combination  of  coefficients  will  be  considered 
admissible  only  if  the  function 

<£*(«, T)  =  I  PV 

satisfies  the  boundary  conditions  at  s=±l;  i.  e.  the  conditions  (b'Kfc*),  and(f')  in  $  2. 
Condition  ( f ' )  (Eq.  (6.  10))  implies 


f  pi :  ^ 

\  kp£,’<>. 


(7.19) 
(7.  20) 

(7.21) 


These  equations  already  imply  the  milder  conditions  (c;),  considering  the  forms  (7.  7), 


(7.9)  of  H  (s). 
n 


Condition  (b'):  From  the  form  (3.  9)  of  j(s)  it  follows  that 

Cn  *■  |  ('  >r  n  (  S+  /  )  dvi(S+  l)  +  0(lsr  H) 

as  s  -•  - 1  ;  hence, 

r(sj  -  OsZ(-')H  (iT,)  +  0(in-ii)  . 


A 

Note  that  constraint  (7.  20)  guarantees  that  no  term  of  the  form  A(s  -  J)  ln(s  -  1)  will 
appear  in  the  expansion  for  s-+l. 

I  am  indebted  to  R.  C.  Ackerberg  for  the  formulation  (7.21). 
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Thus  the  condition  (b')  is  already  satisfied.  Nevertheless  the  asymptotic  equation  (4.2) 
of  Part  I  implies  that  (d/ds)  r^(s)  is  bounded  in  the  neighborhood  of  the  point  s  =  -1. 
Hence  in  the  numerical  computation  we  may  (and  did)  impose  the  additional  constraint 

]  (-Uk  kp2k_  J  =  0.  (7.22) 

k 


29 


8.  NUMERICAL  METHOD 


rA 


We  will  denote 


dt 


and 

t)  p 

=  a~»+  i0/i{c«-’t  e '-H, pi  <-[:(*, P)'1}  hm(i,Zi)hJi,iO-£~rt  <».2> 


These  integrals  converge  if  (7.  19),  (7.  22)  are  satisfied. 

With  the  constraints  (7.  19)  -  (7.  22)  added,  the  variational  problem  implies  by  the 
use  of  the  Lagrange  multipliers  \^,  \^,  X : 

The  coefficients  m  are  listed  in  Table  1. 

nv 


n 

m  . 
nl 

mn2 

mn3 

m 

n4 

=  2k 

1 

0 

k2 

0 

=2k-  1 

0 

k 

0 

(-l)kk 

Table  1 


In  vector  notation  (8.  3)  can  be  written 

B(p)  «•  M  A  -  C  (s.4) 

where  M  is  a  (known)  N  x  4  matrix,  T  a  4-element  column  matrix. 

The  solution  of  (8.  2)  was  obtained  by  Newton's  method.  In  principle,  this  required 
choosing  a  starting  vector  p^  which  satisfies  the  constraints  (7.  19)  through  (7.  22),  and 


30 


to  compute  the  sequences 
relations: 


A  (  jo1"0. 


l*K0) 

P  ,A  , 


(r  =  l,2,3,...)  defined  recursively  by  the  following 


P{t>)  r  M  -  6(  piM)  7  (8.5) 

[>*')  r  0  ,  (8.6) 


where  M  is  the  transpose  of  M.  Equation  (8.  6)  ensures  that  the  constraints  (8.  3) 
continue  to  be  satisfied.  Equations  (8.  5)  and  (8.  6)  represent  a  linear  system  of  N  +  4 
equations,  with  the  unknowns  p^r  +  *-  (n  =  0,  .  .  .  ,  N  -  1)  and  +  (j  =  l,2,3,4).  It 
can  be  shown  that  the  determinant  of  the  system  cannot  vanish  so  that  (8.  5)  -  (8.  6)  always 
has  a  unique  solution. 

The  iteration  for  the  computation  of  the  coefficients  p^  was  programmed  for  the 
IBM  7040  computer  using  this  method.  The  computation  was  started  with  N  =  4,  the 
lowest  number  for  which  the  constraints  can  be  satisfied.  For  N  =  4,  the  coefficients  p 
are  fully  determined  by  the  constraints.  Then  a  number  of  iterations  were  carried  out 
for  N  =  5,  with  the  starting  value  p  =  0;  when  the  function  $(p  ,  .  .  .  ,  p  )  reached  its 

*  O  T 

minimum  for  N  =  5  with  sufficient  approximation,  the  iteration  was  continued  with  N  =  6, 

with  a  starting  value  p,_  =  0,  and  so  on.  Thus,  a  double  iteration  was  used,  consisting  of 

a  sequence  of  iteration  steps  ("major  steps"),  and  each  major  step  was  itself  an  iteration 

for  a  fixed  value  of  N  . 

o 

Practical  details  of  the  computation. 

(1)  In  each  step  of  the  iteration  N(N+  1  )/Z  elements  of  the  matrix  A  and 
N  elements  of  b  have  to  be  evaluated  by  numerical  integration.  In  each  integration  the 
integrand  has  to  be  evaluated  in  M  places,  where  M  is  the  number  of  subdivisions  of 
the  integration  interval  (0,  -j)-  As  a  result,  the  numerical  integration  accounts  for  most 
o.r  the  computer  time  needed  to  perform  the  Newton  iteration.  We  used  in  the  integrations 
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v 


-y 


7 

Gauss  mechanical  quadrature,  96-point  formula.  The  integration  subroutine  yielded 

5  or  more  significant  digits  accuracy,  when  the  integrands  were  bounded,  reasonably 

smooth  functions.  To  test  the  integration  routine,  a  ,  b  were  evaluated  with  numerical 

mn  n 

integration,  and  the  results  were  compared  with  the  closed  formulas  derived.  (Appendix  3) 
(2)  The  error  of  each  step  of  Newton's  iteration  was  measured  by 

C  =  ,wax  I  B  ( r> )  I  /  max  |  5.  I 

k  *  '  k  k 

where  is  the  biggest  partieil  sum  (in  absolute  value)  in  the  computation  of  B^.  The 
iteration  in  a  major  step  was  considered  finished  when  e  did  not  exceed  2x  10  It 
was  somewhat  surprising  that  with  the  exception  of  the  first  three  major  steps,  this  limit 
was  reached  in  a  single  iteration.  At  the  end  of  each  major  step,  we  computed  the  new 
approximation  to  and  the  error  6(t)  committed  in  the  boundary  condition  (5.  1),  in 

the  96  subdivision  points  of  integration,  i.  e.  , 

S(t)  =  iknh  C\  Uv,k  R(  l, it)  f  (l,  It) 

-  [Ik  2  I  >„  KJ 


The  square  norm  of  6(t). 


i$h 


was  used  to  judge  the  goodness  of  the  answer  in  each  major  step.  As  an  additional  check 
the  total  angle  of  deflection  of  the  jet  was  computed  after  each  major  step  using  the  for¬ 


mula 


- 

0mflx  “  )  QCY,0) 

o 

--I”7'1  {<***  E(t,p)~  E(t,p) }  4£-  - 

J  0 


(The  theoretical  limit  value  of  &max  after  infinitely  many  iterations  is  of  course  90  ). 


The  values  of  $ ( p) ,  6  ,  P  obtained  in  the  course  of  the  iteration  are  listed  in 

Table  2.  As  is  evident  from  this  table,  ||6||  decreased  steadily  during  the  first  11 

-4 

major  steps,  and  reached  a  minimum  of  4.  38  x  10  ,  well  within  the  tolerance  of  engi¬ 

neering  computations.  After  this  the  goodness  of  the  approximation,  as  judged  from  the 
values  of  ||6||  ,  deteriorated  first  slowly,  then  faster,  and  in  general  changed  in  an 
erratic  fashion.  The  reason  for  this  peculiar  trend  is  probably  that  the  integration  routine 
loses  accuracy  for  high  harmonics  sinnT,  cos  nT,  while  the  theoretical  gain  in  accuracy 
of  these  terms  is  becoming  smaller.  Also,  the  SHARE  routine  used  for  the  solution  of 
systems  of  linear  equations  loses  accuracy  for  a  large  number  of  equations.  Another 
measure  of  the  progress  of  the  computation  is  the  decrease  of  $(p).  As  seen  from  Table 
2,  f(p)  converges  extremely  rapidly  to  its  minimum,  in  fact  much  more  rapidly  than  6. 
(The  best  value  of  $(p)  was  0.474498,  whereas  the  first  (!)  major  step  resulted  in 
$(p)  -  0.  480.  .  .  )  The  coefficients  A.  .(p)  also  converge  rapidly;  they  hardly  change  from 
one  step  to  the  next,  once  computed.  However,  the  coefficients  p^  do  not  show  any  trend 
of  convergence,  a  fact  explained  by  the  infinitely  many  possibilities  to  represent  a  single 

function  Q  (1.  r)  by  a  series  p  h  (1,t).  The  "best"  choice  of  the  coefficients  p 
o  rn  n  n 

(n  =  0,  1,  .  .  .  ,  14)  is  listed  in  Table  3. 


The  iteration  described  above  took  5^  minutes  on  the  IBM  7040  computer. 

(3)  In  the  final  phase  of  the  computation  the  coefficients  p^  were  used  to 
obtain  the  functions  X (cp,  l|f)  and  Y(cp,i|f).  This  was  based  on  the  formula 


z-  cj  ~  7  ►  c  [w-2(-w)'~] 

Jo  L  L  j 


(8.  7) 


where 


C  =  e 


-  r  A  -  1 ) 
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to 


‘■ft 


ok 


Equation  (8.  7)  is  a  direct  consequence  of  Part  I,  Eq.  (4.  5).  This  form  of  the  equation 
has  the  advantage  that  the  integrand  remains  bounded  at  the  origin.  The  integration  was 
carried  out  by  a  subroutine  to  solve  systems  of  ordinary  differential  equations,  based  on 
the  Adams-Bashford  method  (4^  order),  written  by  Kenneth  Plotkin. 


N 

V 

$(p) 

INI 

0  /degr. 

max/ 

T^y  sec. 

5 

5 

0.48025144 

3.  807  x  10"1 

66.  58 

.  68 

6 

2 

0.  47479781 

1.  671  x  10_1 

93.  664 

.  90 

7 

2 

0.47450213 

2.  347  x  10‘2 

87. 932 

1.  17 

8 

1 

0.  47449958 

3.  823  x  10‘3 

88. 223 

1.  48 

9 

1 

0.  47449929 

7.  608  x  10"3 

88.  6648 

1.  83 

10 

1 

0.  47449854 

1.631  x  10‘3 

89.  5812 

2.  23  > 

11 

1 

.  0.  47449850 

4.  156  x  10"3 

89. 6775 

2.  68 

12 

1 

0.  47449845 

1.  017  x  10‘3 

89.  7799 

3.  18 

13 

1 

0.  47449843 

4.  63  x  10-4 

89.992 

3.  73 

14 

1 

0.  47449842 

4.  62  x  10-4 

89. 968 

4.  32 

15 

1 

0.  47449841 

i 

o 

X 

00 

rO 

rf 

89.  974 

4.  97 

1 

16 

1 

0.  47449841 

1.  676  x  10’3 

90.  027 

5.  65 

17 

1 

.  0.  47449841 

1.  587  x  10‘3 

89.  979 

6.  37 

18 

1 

0.  47449839 

i 

o 

X 

oo 

oc 

u . 

89. 921 

7.  17 

19 

1 

0.  47449838 

1.  640  x  10"3 

89.  889 

8.  00 

20 

1 

0.  47449838 

2. 257  x  10’3 

89. 899 

8.  88 

Table  2. 


N  number  of  unknown  coefficients  p 

rn 

v  number  of  iterations /major  step 

<{>(p)  value  of  variational  integral 

1 1 6  ||  norm  of  error  in  boundary  condition 
9  total  angle  of  deflection  of  jet 

Tj  time  required  for  individual  iteration 
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n 

/ 

Pn 

n 

Pn 

0 

. 520853172  x  10_1 

1 

433968985  x  10‘2 

106932819  x  10'2 

2 

261440941  x  10_1 

158038496  x  10'1 

3 

225389441  x  10'1 

210494106  x  10'2 

m 

410373311  x  10-2 

. 102374660  x  10"1 

5 

. 461629598  x  10'2 

. 168663154  x  10‘2 

6 

. 676816046  x  10-3 

155102682  x  10'2 

n 

251967831  x  10-3 

149858572  x  10‘3 

Table  3. 


'Pin  P^Pin-l-Pn* 
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9.  DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 


Plots  of  the  bounding  streamline  which  were  determined  numerically  using  a 

variational  technique  are  displayed  in  Figs.  3  and  5,  along  with  Taylor's  theoretical 

1  3 

results  and  some  experimental  data  from  Vizel  and  Mostinskii  .  Some  of  the  stream¬ 
lines  and  equipotential  lines  are  shown  in  Fig.  6  in  the  neighborhood  of  the  jet  opening. 
The  numerical  results  indicate  a  deeper  jet  penetration  than  had  been  found  by  either 
Taylor  or  experiment.  Taylor  mentions  the  difficulty  of  verifying  theoretical  results 
experimentally  due  to  the  viscous  spreading  of  the  jet  which  would  fill  a  wedge  of  nearly 
40  degrees.  With  this  in  mind,  the  discrepancy  between  the  theoretical  curve  and  the 
experimental  data  in  Fig.  5  is  probably  not  so  great. 

Finally,  using  Bernoulli's  principle  in  the  external  flow,  the  coefficient  of  the 
pressure  can  be  written 


c 


f 


p  ~  )°g- 


I  - 


(9.1) 


Numerical  values  for  C  were  computed  versus  distance  along  the  plate  upstream  and 

P 

are  displayed  in  Fig.  7.  It  is  readily  seen  that  the  pressure  decreases  slowly  from 
the  stagnation  pressure  at  the  jet  exit  to  the  free  stream  value  many  slot  widths 
upstream. 
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APPENDIX  1 


Proof  of  Equation  (5.  4): 


+  >o 


To  show  (5.  4)  we  transform  A^  Lh,  k]  by  Green's  formula,  taking  into  account 

€ 

that  h(cp,  i|i)  is  harmonic.  Thus, 

D*[K,fcJ=-  lIL  VhVKd<fdi  -  if  ^<<5, 

t  C*  (.T  * 

tfP;, ■  P h  < 1  k  H s 
5  I  {  C  '  C) 


and 


(Al.  1) 


where 


f  *  S.  -I L<1  xkd*>  ) 

c  \  L  (■ 

=  d,-(T'»  V1  ds  ' 


Let  now  e-0.  Then,  using  property  (b),  a  simple  computation  shows  that  X-0, 
and  since  k  is  bounded, 


/ 


k  cl 


C  . 


r 

Letting  X-oo,  from  property  (c/  follows  that  I  kds  -  0.  Substituting  these  limits  into 


(Al.  1),  the  statement  (5.  4)  follows. 


38 


APPENDIX  2 


Proof  of  the  identity  (5.  6): 


AL«J=  l  , 

A  '■A 

X  r  i  #  * 

A  LuJ  can  be  written  by  the  substitution  u  =  u  t  -  (n  r,  using  polar  coordinates 


(  r,  9),  in  the  form 


| 

•(  i 

I  <L*\  *. 

v 2  . 
)  4  l 

LaJ 

\ 

V  l 

.  \  c  -t  n-t 

J  4  M 

■  se  '  J 

IT  >2 

f  j  t  % i«0  j  = 

^  7T 


2  7T  (a\  f 

r* 

-  [{  (va');toihletf  ]  uft(A,e)<*e 

J  **  A  '  m 


+ - A;  -  f  f  <■(.'(  f  ,  6  )  (J 

2 n  2v[€  L  J._ 


e 


-+■  ^Vi  A  . 

JT 


Let  now  €-0.  Taking  into  account  that  u  is  bounded  in  the  neighborhood  of  the  origin, 
we  obtain 


4a"[«J-  40'V'J  r  T  j.t  ^  %  A  . 

w  from  the  identity 

Af«  ]  A*[«J  *  DAC«]  - 


The  proposition  follows  now  from  the  identity 

A 
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APPENDIX  3 


Proof  of  the  convexity  of  the  functional  F[u]. 


We  first  show  that  the  (ordinary)  Dirichlet- integral  is  a  convex  functional.  Let 

S  denote  an  arbitrary  domain,  D  [u]  the  Dirichlet-integral  of  the  function  u  in  S  and 

d[u]  =  (D  [u])^,  the  "Dirichlet- norm"  of  u.  If  d[u]  <  oo,  d[v]<  cd,  then  u,  v  satisfy 
8 

g 

the  following  the  triangle  inequality 

cil><-vj£r([u3  '  <A3-‘> 

Let  k  by  any  number,  0  <  k  <  1.  Then,  since  d[u]  is  homogeneous  first  order 
#U)  =  ku  f-(  I-  <.)v]  <  k  d[«-]  t(l~k)o([v]  ,  (A3.2) 

i.  e.  ,  d[u]  is  convex  (but  not  strictly,  since  equality  is  possible  in  (A3.  2)  ).  Moreover, 

since  u,  v  are  arbitrary,  $(k)  is  a  convex,  non-negative  function.  Actually,  $(k)  can 

vanish  at  most  for  a  single  value  of  k,  for  which  ku+(l-k)v  =  0  almost  everywhere 

2 

(possible  only  if  u  and  v  are  linearly  dependent).  This  implies  that  [$(k)]  is  strictly 
convex;  consequently,  if  u  ^  v, 

(A3.  3) 


i.  e. 


Ds[kut  K  DsL«jKl-k)DsC^  ' 

Apply  now  (A3.  3)  to  S=G^,  and  the  admissible  functions  u,  v,  u  ^  v.  and  also  to 

\  *  :'c  r 

S  =  G  ,  u  ,v  .  Then  from  (5.  6)  follows  the  convexity  of  ALuJ. 


The  convexity  of 

R[u]  -  \  [  CtriA  UH.O)-  t 

^  0 

is  an  immediate  consequence  of  the  convexity  of  the  function  cosh  u.  Therefore 
F[u]  =  a[u]  +  B[u]  is  also  convex. 
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APPENDIX  4 


Evaluation  of  the  Matrices  (a  ),  (b  ). 

mn  n 


a  ,  ,  =  0  i  i  rn  f 

2  m  ,2  n 

_  7T 


M  , 


a 


2  n  ?  2  n 


2  m ,  2 


2  m 


»>'-*  m 


(X^  =  0  if  k  is  even,  =1,  if  k  is  odd) 


Q  -  _  £  m  n  / 

2m  - 1 ,2n  -i  =  ^  ‘ 


-cc>n)  >f  n  +  n. 


a 


2n 


,  -  o 

~  7T  rv  ' 


Here 


1 ♦  i  f 

t0<2A  *  I  ~  1 

=  I*  i** 

Zk  *  i  ~  l  *  } 


.  .  .  *  > 
2  k  '  I 


r±-L 


2k  I  ?■  2  k  *  I 


*  (2  k  «/ 


f 


*  *T"  .  ' 


(2k  -  I )*~  l  (2  k*  0^ 

The  matrix  a  is  symmetric,  i.  e.  ,  a  =  a  For  both  even  and  for  odd  subscripts 

mn  nm 

only  non-negative  values  may  be  substituted. 


The  matrix  b  : 

*>c  =  f  '  =  0  'f  k*C  ■ 

^Ak  *■  1  2k  k  i  *  %Xr  5  ^ 
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FIG.  la.  REGION  OF  FLOW  IN  PHYSICAL  PLANE 
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FIG.  3.  THE  BOUNDING  STREAMLINE 
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FIG.  6.  STREAMLINES  AND  EQ UIPOTENTIAL  LINES 


